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class (M ~ 50M o ), in a phase of moderately super-Eddington accretion. The formation of such a massive BH via normal stellar 
evolution may have been favoured by the very metal-poor environment of NGC 5408. 



Abstract. We have studied an ultra-luminous X-ray source (ULX) in the dwarf galaxy NGC 5408 with a series of XMM- 
Newton observations, between 2001 July and 2003 January. We find that its X-ray spectrum is best fitted with a power law of 
1 photon index T ss 2.6-2.9 and a thermal component with blackbody temperature kT^ ss 0.12-0. 14 keV. These spectral features, 

and the inferred luminosity » 10 40 erg s -1 in the 0.3-12 keV band, are typical of bright ULXs in nearby dwarf galaxies. The 
blackbody plus power-law model is a significantly better fit than either a simple power law or a broken power law (although 
^ ' the latter model is also acceptable at some epochs). Doppler-boosted emission from a relativistic jet is not required, although 

we cannot rule out this scenario. Our preliminary timing analysis shows flaring behaviour which we interpret as variability in 
G\ . the power-law component, on timescales of ~ 10" s. The hard component is suppressed during the dips, while the soft thermal 

component is consistent with being constant. The power density spectrum is flat at low frequencies, has a break at Vb ~ 2.5 
mHz, and has a slope w -1 at higher frequencies. A comparison with the power spectra of Cyg X-l and of a sample of other 
^-r 1 BH candidates and AGN suggests a mass of ~ 10 2 M G . It is also possible that the BH is at the upper end of the stellar-mass 

o 
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a. 

c/3 1. Introduction & Rees 2001), or runaway stellar mergers in young clusters 

(Portegies Zwart et al. 2004). 

> ! Ultra-luminous X-ray sources (ULXs) are point-like sources, A ULX in the starburst dwarf irregular galaxy NGC 5408 

"£j ; not including galactic nuclei and young supernova remnants, could help discriminate between these two alternatives. Radio 

■ with apparent luminositites higher than the Eddington limit observations of this galaxy (Stevens et al. 2002) showed four 

C3 | for a stellar-mass accreting black hole (BH), ie, with L x £ ma i n emission regions, mostly coincident with massive young 

3 x 10 39 erg s" 1 (Fabbiano 1992; Colbert & Mushotzky 1999; star c i uste rs and Ha emission. The X-ray emission is domi- 

Roberts & Warwick 2000). The true nature of these objects nated by a single point source , with a luminosity > 10 40 erg 

is still open to debate. The fundamental issue is whether the s -i (Fabian & Ward 1993; Kaaret et al. 2003), located outside 

emission is isotropic or beamed along our line-of-sight; in the HII reg ions, but apparently coincident with a weak, steep- 

the latter case, the total luminosity would of course be lower spe ctrum radio source. On the basis of the X-ray, radio, and op- 

than the isotropic value. A possible scenario for geometri- tical fluxeSi Kaaret et al (2 003) concluded that the most likely 

cal beaming is a super-Eddington mass inflow during phases explanation for the ULX was beamed relativistic jet emission 

of thermal-timescale mass transfer (King 2002). Relativistic from a st ellar-mass BH. However, they could not rule out an 

beaming, associated with Doppler boosting, would instead be intermediate-mass BH model. In this paper we present the first 

the effect of our looking into the jet of a microquasar (Kording results of our XMM-Newton study of the ULX in NGC 5408, 

et al. 2002). Alternatively, if the emission is isotropic and the providing new insights into the nature of this source. 
Eddington limit is not violated, ULXs must be fuelled by ac- 
cretion onto intermediate-mass BHs, with masses ~ 1O 2 M . 
Possible mechanisms of formation of such massive remnants 
include the collapse of massive population III stars (Madau 



2. Observations and data analysis 



NGC 5408 was observed with the European Photon Imaging 

Camera (EPIC) and the Reflection Grating Spectrometer 

Send offprint requests to: R. Soria (RGS) on board XMM-Newton in 2001 July-August, 2002 July 
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and 2003 January (Table 1). The Optical Monitor was blocked 
because of a bright foreground star in the field. In this paper, we 
present preliminary results from the EPIC data; more detailed 
results are left to further work. 

The thin-filter, full-frame mode was used for the MOS and 
pn cameras. We processed the Observation Data Files with 
standard tasks in Version 5.4 of the Science Analysis System 
(SAS). After inspecting the background flux of each expo- 
sure, we rejected various intervals with background flares in 
Rev. 301 and Rev. 574. A high, flaring background affected the 
whole exposures in Rev. 483: the average background count 
rate was « 20% of the ULX count rate in the source extraction 
region. We retained that exposure (the only one available for 
2002), with the caveat that the timing analyses may be affected 
by larger errors. For all other exposures, the background rate 
in the selected good-time-intervals is typically ~ 1-2% of the 
source count rate. Technical problems during Rev. 313 meant 
that the pn was not activated, while the usable part of the MOS 
exposures lasted only 6.9 ks despite a nominal exposure time 
of 7.5 ks. We filtered the event files, selecting only the best- 
calibrated events (pattern < 12 for the MOSs, pattern < 4 for 
the pn), and rejecting flagged events. We also checked that the 
source was not piled-up, with the SAS task epatplot. 

For each exposure, we used a 30"-radius source extrac- 
tion region, with the background extracted from other re- 
gions in the same CCD. We built appropriate response func- 
tions with the SAS tasks rmfgen and arfgen. We fitted the 
background-subtracted spectra with standard models in XSPEC 
v.11.3 (Arnaud 1996); owing to the uncertainties in the EPIC 
responses at low energies, we used only the 0.3-12 keV range. 
Firstly, we analysed the individual pn and MOS spectra from 
each observation. After ascertaining that they were consistent 
with each other, we coadded them with the method described 
in Page et al. (2003), to increase the signal-to-noise ratio. Thus, 
we obtained one combined EPIC spectrum of the ULX from 
each epoch. With the same method, we also built an aver- 
age EPIC spectrum of all the observations from 2001. For the 
timing analysis presented here, we selected all the good-time- 
intervals when both the MOS and pn were operating (except 
for the lightcurve in Rev. 313, for which only the MOS is 
available), and we retained the events in the 0.2-12 keV range. 
We used standard tasks in FTOOLS (Blackburn 1995) to obtain 
background-subtracted lightcurves and power-density spectra. 



3. Main results 

3.1. Spectral analysis 

Based on a short (4.7 ks) Chandra observation, it was sug- 
gested (Kaaret et al. 2003) that the X-ray spectrum of this ULX 
could be fitted equally well by a broken power law (break en- 
ergy « 0.65 keV), or by an absorbed blackbody plus power- 
law model (kT\,b « 0.11 keV; photon index F « 2.8). The 
broken power-law model would support the interpretation of 
this source as beamed emission via inverse-Compton scatter- 
ing of optical/UV photons by a relativistic jet (Georganopoulos 
et al. 2001; Kaaret et al. 2003). 



Table 1. XMM-Newton/EPIC observation log for NGC 5408. 



Date 


Obs ID 


Start time Stop time 


GTI 


2001 Jul 31 


0112290501 


pn. LZ.^y .t+j 


i a- 1 ^-n^ 


3.8 ks 




(Rev. 301) 


MOS: 12:12:08 


14:17:21 


4.1 ks 


2001 Aus 08 


0112290601 


pn: 10:22:14 


11:45:34 


4.5 ks 




(Rev. 305) 


MOS: 09:43:01 


11:48:13 


7.4 ks 


2001 Aug 24 


0112290701 


pn: - 




- 




(Rev. 313) 


MOS: 21:58:47 


00:05:39 


6.9 ks 


2002 Jul 29 


0112291001 


pn: 08:23:12 


10:08:22 


4.9 ks 




(Rev. 483) 


MOS: 08:01:02 


10:11:28 


7.7 ks 


2003 Jan 28 


0112291201 


pn: 00:06:06 


01:29:57 


3.0 ks 




(Rev. 574) 


MOS: 23:43:56 


01:18:14 


3.3 ks 



In our XMM-Newton study, we confirm that a simple ab- 
sorbed power-law model is ruled out (x\ ^1.4 for all spectra). 
We then consider a broken power-law model; we impose the 
physical constraint that the photon index below the break must 
be lower than the index above the break. Although at some 
epochs this model provides a good fit, with parameters simi- 
lar to those of Kaaret et al. (2003), at other epochs (e.g. for 
Rev. 305, 2001 Aug 08) it leads to large, systematic residuals 
(Tables 2, 3). Instead, a blackbody/disk-blackbody plus power- 
law model provides a consistently good fit over all the obser- 
vations. In this model, the power-law index T a 2.6-2.9 is sim- 
ilar to that found in Galactic BH candidates in the high/soft 
and very high states (Belloni 2001). The low-temperature ther- 
mal component is similar to the soft emission found in bright 
ULXs, such as those in NGC 1313 (Miller et al. 2003), Ho II 
(Dewangan et al. 2004) and NGC 4559 (Cropper et al. 2004). 
We cannot discriminate from our observations between a 
pure blackbody and a multicolor disk-blackbody spectrum. 
Comptonisation models such as brae (Shrader & Titarchuk 
1999) also provide a good fit (Tables 2, 3). The thermal com- 
ponent is significantly detected at all epochs, with similar tem- 
peratures (kT bh » 0.12-0.14 keV). We find long-term flux and 
spectral variations over the 18-month span of our observations 
(Fig. 1 and Tables 2, 3), but no real state transitions. At an 
assumed distance of 4.8 Mpc (Karachentsev et al. 2002), the 
emitted luminosity in the 0.3-12 keV band is « 10 40 erg s _1 . 

In the spectral models described above, we have assumed a 
solar metal abundance for the intrinsic absorber (wabs model in 
XSPEC). In fact, the metal abundance of the starburst region in 
NGC 5408 is known to be lower, Z * 1.2 x 10~ 3 ~ 1/14Z 
(Stewart et al. 1982). It is therefore likely that the intrinsic 
absorber in front of the ULX has also low metal abundance. 
We have repeated the fits using a tbvarabs model (Wilms et 
al. 2000) for the intrinsic absorber. Firstly, we left the metal 
abundance as a free parameter: we found that its value is not 
well constrained, and we did not obtain significant^ 2 improve- 
ments for any of the models. This is because the intrinsic col- 
umn density is low, hence the fitted models are not much af- 
fected by changes in its metal abundance. We then fixed the 
metal abundance at Z = 0.07Z G , and repeated the fits for 
the broken-power-law (bknpo) and blackbody-plus-power-law 



R. Soria, C. Motch, A. M. Read and I. R. Stevens: Flares from the ULX in NGC 5408 



3 



(bb+po) models (Tables 4, 5). Even with a lower metal abun- 
dance, including a low-temperature thermal component signif- 
icantly improves the fits, compared with a simple power law. 
We also find once again that bb+po models give a better xl 
than bknpo models, for 4 out of the 5 spectral observations, 
although we cannot rule out the broken-power-law model alto- 
gether. 

We tried to find a quantitative way to estimate whether the 
bb+po model is "significantly" better than the bknpo model, as 
suggested by the differences in xl an d null hypothesis proba- 
bilities, or whether the differences between the two can be at- 
tributed to chance. Two models can be directly compared with 
the F-test (Hamilton 1965) in the special case when the param- 
eters of one model are a subset of the parameters of the other 
model. In a more general case, such as the one we are dealing 
with, a quantitative test for comparing the predictions of two 
models with the observed values in each data bin was first dis- 
cussed in Williams & Kloot (1953), and refined in Himmelblau 
(1970) and Prince (1982). We summarize the results of our 
comparison in Table 6: in 4 of the 5 observations, there is a 
^ 90% probability that the bb+po model is a significantly better 
fit than the bknpo model. Only for the last observation (2003 
Jan) are the two models statistically indistinguishable. 

A bb+po model is also a significantly better fit than the 
bknpo model for the average of the 3 observations from 2001 
(Tables 2-5). This is not in itself a proof against the bknpo 
model: it may simply result from the individual spectra being 
broken power-laws but with different parameters. 

In conclusion, we argue that a bb+po model, or any other 
similar model with a power-law and a soft thermal component 
(eg, a Comptonisation model such as bmc), are significantly 
better fits to the XMM-Newton data than either a simple or a 
broken power law. This does not rule out a relativistic jet in- 
terpretation for this ULX: it was shown by Georganopoulos et 
al. (2001) that a broken power-law spectrum is only an approx- 
imation, and that the true spectrum would in fact be curved 
around the break energy (in our case, « 0.5-0.6 keV). Such a 
spectral model would be practically indistinguishable from a 
power-law plus thermal component: observations over a larger 
energy range are needed to test this prediction. 

3.2. Timing analysis 

Short-term variability, significant to > 99% according to the 
Kolmogorov-Smirnoff and x 2 tests, is detected at all epochs, 
except during Rev. 313 (which has a lower signal-to-noise ra- 
tio, since the pn camera was not operating). The source exhibits 
a flaring behaviour, with changes in flux by a factor of ~ 2 
over ~ 100 s (Fig. 3). The flares are most evident in the 2003 
observation, with a Kolmogorov-Smirnoff probability of con- 
stancy < 10~ 3 . Here we focus on two important results from 
this epoch; a detailed timing study of all the other observations 
will be presented elsewhere. 

Comparing the background-subtracted lightcurves in a soft 
(0.2-1.5 keV) and hard (1.5-12 keV) band, we found that the 
hard flux drops to a value consistent with zero during the dips 
(Fig. 4, top); in the soft band, a "baseline", persistent corn- 



Table 2. Best-fit parameters for the combined EPIC spectra of 
the ULX during the three observations in 2001 (Rev 301: Jul 
31; Rev. 305: Aug 08; Rev. 313: Aug 24). The quoted errors 
are the 90% confidence limit for one parameter (Ax 2 = 2.7). 
We assumed a Galactic column density n H ,Gai = 5.7 x 10 20 
cm -2 (Dickey & Lockman 1990). Fluxes and luminosities are 
the unabsorbed values. 



Parameter 


Value Jul 31 


Value Aug 08 


Value Aug 24 


model: wabs Ga i x wabs x bknpo 


r, 
1 1 

£ br (keV) 

r 2 

^ bp (io- 3 ) 


2 +03 
98+ 81 
55 +004 

-0.04 

Q TO+0.14 

-0.13 

3 5 + " 


2 +03 

"■^-0.2 
i 79+0.47 
1 -' z '-0.58 

63 +0 05 
u - UJ -o.io 

■X ^4+0.13 
-0.12 

1 7 +2 ' 7 


< 0.3 

65 +L2S 

"• UJ -1.44 

53+ 013 

-0.05 
•2 i c+0.15 
J -0.07 

3.3!" 


xl 

/0.3-12 (b) 
A).3-12 <C) 


1.16(164.2/141) 

3.63-j 
1.0 


1.45 (229.0/158) 
3 5 +0 ' 2 
1.0 


0.99(132.2/133) 
3.2!« 4 
0.9 




model: wabs Ga i 


x wabs x (bb + po) 


«H <a) 

kT bh (keV) 

K bb (io- 5 ) 
r 

K po (10- 4 ) 


4 +03 

u -^-0.2 

14 +001 

u - -0.02 

2 5 +a9 
9 q 1+020 

^•^ -0.19 

5 7 +L2 
J -'-l.l 


3 +0 ' 2 

15 +001 
u - 1J -o.oi 

9 Q+0.5 

2 65 +018 
4 6 +a9 


0.13t°> 

3 2 +2 - 3 
J - z -i.o 

O CC+0-16 
/ - 88 -0.18 

6 7 +l 3 


xl 

f (b) 
/0.3-12 

Lfl.3-12 <C) 


1.00(141.4/141) 
1.3 


1.09(172.0/158) 

t a +0.2 

1.1 


0.86(113.9/133) 
5 5 +0 - 3 
1.5 


model: wabs Ga i x wabs x (diskbb + po) 


kT m (keV) 

^dbb 

r 

£po (io- 4 ) 


n £+0.2 
17 +002 

1 ' -0.02 
977+390 
Z "-143 
9 85+0.19 
z "° -0.20 

5 3 +1 ' 3 


5 +0 ' 2 

" -0.2 

19 +001 

U - I: *-0.02 
1?4 +178 

1 '^-48 
9 cc+0.24 

-0.17 

4 1 + 1 ' 2 
H --0.7 


1 l +05 
15 +003 

J -0.02 
ccc+1580 
JJJ -456 

2-86!^ 
6-6^ 


xl 

/0.3-12 (b) 
A).3-12 <C) 


0.99(139.6/141) 

5.2!?: 6 
1.4 


1.05(166.2/158) 
4 5 +(u 
1.2 


0.86(114.4/133) 
1.8 




model: wabs, 


3a i x wabs x bmc 




RH W 

kT bb (keV) 

r 

log(A) 

^bmc (IO" 5 ) 


< 0.4 
13 +001 

-0.01 

2.94*"' 
-0.18!™ 

3 9 +L3 


< 0.2 

2.68^ 
-0 30 +0 12 

u - JU -0 .12 
t 4+O.6 


5 +05 
12 +0 - 02 
2 90 +0 - 22 
-0 24 +019 

u - z +_0.16 

5 3 +24 


xl 

f (b) 
/0.3-12 

A).3-12 <C) 


1.02(144.2/141) 
3 4 +al 

-'•^-0.8 

0.9 


1.13(178.4/158) 

o 1+0.4 
-'• 1 -0.8 

0.9 


0.86(114.8/133) 
44+1.1 

^•^-1.7 

1.2 



b in units of 10 12 erg cm 2 s 1 
c in units of 10 40 erg s~' 
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Table 3. Best-fit parameters for the combined EPIC spectra of 
the ULX in 2001 (Jul 24-Aug 25), 2002 (Jul 29) and 2003 (Jan 
28). The quoted errors are the 90% confidence limit for one 
parameter (A^ 2 = 2.7). We assumed a Galactic column density 
«H,Gd = 5.7 x 10 20 cnr 2 (Dickey & Lockman 1990). Fluxes 
and luminosities are the unabsorbed values. 



Table 4. As in Table 2, but with sub-solar abundances (Z 
O.O7Z ) for the intrinsic absorber. 



Parameter 


Value in 2001 


Value in 2002 


Value in 2003 




model: wabs Ga i 


x wabs x bknpo 




r, 

£ br (keV) 

r 2 

K bp do- 3 ) 


r° 1 

76 +0 16 

-0.48 

56 +0 05 

" -0.03 
Q 9 1 +0.03 
-0.02 

■x 4+0.2 

J - -0.2 


< 0.2 

°- 58 -U5 

60+ 005 

"• u -0.06 

2 89 +o u 
2 +3 3 

■^""-0.6 


< 0.2 

1 1Q+0.68 
1 ' 1 '-1.75 

55 +0 05 

J -0.03 

2.87*»* 
6-2^ 


xl 

/o.3-i2 (b) 

1*3-12 <C) 


1.79 (245.4/137) 

i 9+O.2 
J - -0.2 

0.9 


1.09(99.2/91) 

9 i+O.l 

^-o.i 
0.6 


0.97 (83.2/86) 
2 +0 ' 

^•"-0.1 

0.6 


model: wabs Ga i x wabs x (bb + po) 


m <"> 

kT bh (keV) 

Kbb (io- 5 ) 
r 

K vo (IO" 4 ) 


"'^-0.2 

14 +001 
2 5 +05 

-0.3 
9 7Q+0.11 

■'•"-o.io 

c 7+O.6 
J -'-0.6 


9 +0 1 
13 +0 - 01 

-0.02 
9 c+0.4 
^• J -0.5 
9 co +0.06 
/ --> 8 -0.06 

5 +0 - 2 

J - -0.3 


i 1+0.6 
1 - J -0.5 

12 +002 
"■^-0.02 

3.1 +3 - 5 

2.6811? 
4.9_J;° 


f (b) 
70.3-12 

Lfl.3-12 <C) 


1.10(150.9/137) 
4 5 +01 
1.2 


0.98 (88.9/91) 
4 +01 

1.1 


0.99 (85.5/86) 

4 -x+a.2 

1.2 



Parameter 


Value Jul 31 


Value Aug 08 Value Aug 24 


model: wabs Ga i x tbvarabs x bknpo 


n (») 
£ br (keV) 

r 2 

^ bp (io- 3 ) 


U -°-0.4 

1 61 +1 30 
53 +0 06 

-0.06 

3 40+°- 14 

J - H "-0.13 
9 7+5.8 
Z ''-1.5 


n q+ - 4 n ?+ - 7 

2.12:;- 68 i-02 + ;^ 

51 +0 - 35 57+ 013 

J -0.20 -0.13 

T, 00+0.10 O 1 q+0.15 

J - JO -0.11 -'• IO -0.07 

2.0!: 2.8« 


f (b) 
70.3-12 

£fl.3-12 <C) 


1.15(162.3/141) 

4 o+ - 8 

^■"-0.7 
1.1 


1.43(225.3/158) 0.99(132.0/133) 

4 4+10 o 4 +1.0 
H - -1.1 -0.2 

1.2 0.9 




111UUC1. WaUsQyl 




«B W 

££ b b (keV) 

^bb (io- 5 ) 
r 

K po (io- 4 ) 


5 +02 

—0.2 

14 +001 
u - iH -o.oi 

1 O+0.4 
i -'-0.5 

9 ni+0.17 
z " 7J -0.23 
c 7+0.9 
—1.2 


3 +0 2 7 +a2 
15 +001 13 +001 

"• 1J -0.01 "- 1J -0.02 

2 +02 1 9 +1 3 

-0.3 ' -0.9 
9 £]+0.18 9 oc+0.16 
-0.18 ^- OJ -0.18 
4 4+0.9 fi 4+1.0 


/0.3-!2 (b) 
£■0.3-12 <C) 


1.00(141.6/141) 

4 1+0.5 

1.1 


1.10(173.2/158) 0.87(115.2/133) 

•3 5+0.3 4 4+O.5 
J - -0.2 H - -0.5 

1.0 1.2 



a in units of 10 21 cm -2 

b in units of IO -12 erg cm -2 : 

c in units of 10 40 erg s -1 



model: wabs Ga i x wabs x (diskbb + po) 



«H <a) 

kTm (keV) 

^dbb 

r 

Kpo (IO" 4 ) 


7 +a2 
17 +001 

'-0.01 
97^ + 165 

9 70+O.H 
z - /J -0.11 

5 3 +a7 


1 3 +05 

15 +001 
u - IJ -o.oi 

471+952 

-289 

9 C7+O.I9 
z.j;_ on 

5 +09 


1 5 +0 - 7 
15 +002 

"■ iJ -0.02 

670: 23 °° 

2.66^1? 
4 8 +L1 


xl 

/0.3-12 (t) 
£-0.3-12 ° 


1.05(144.5/137) 

5 2 +a2 
J - -0.8 

1.4 


1.02 (93.1/91) 
5 +2 
1.4 


1.02 (87.4/86) 
c 9+4.0 

1.4 




model: wabs Ga i 


x wabs x bmc 




kT bh (keV) 

r 

log(A) 

^bmc (IO" 5 ) 


1 +a8 
0-13^1 

-0.20!«i° 

J.8_ 02 


0.6*}| 
0.13t° 2 
2 60 +0 - 07 

^■•""-0.07 

-0 24 +0 - 04 

"•■^-0.04 

3 6 +0 - 3 


9 +0 - 4 

"•^-0.4 

12 +002 

u - lx -0.02 

2 67 +a20 

•^•"'-0.17 

-0 30 +0 20 
4 +35 

^•"-2.0 


f (b) 
70.3-12 

£0.3-12 <C) 


1.17(160.5/137) 

3 5 +a3 
1.0 


0.95 (86.4/91) 

0+0.4 

J - J -1.4 

0.9 


0.97 (83.5/86) 
4+0.4 

0.9 



3 in units of 10 21 cm -2 

b in units of 10 -12 erg cm" 2 s 

in units of 10 40 erg s" 1 



ponent is always present in addition to the flaring component 
(Fig. 4, middle). The rms fractional variability is « 70% in the 
hard band, and « 30% in the soft band. There is no time lag be- 
tween the flares in the two bands, suggesting that we are seeing 
the same variable physical component in both. 

The flaring behaviour in the two bands has a simple in- 
terpretation in the framework of the bb+po spectral model 
(Sect. 3.1). In the hard band, the variability is necessarily due 
to the power-law component: there is no contribution from the 
thermal component in this band. A change in the power-law 
component would also affect the soft band, where the power- 
law and thermal fluxes are of the same order of magnitude. 
Hence, we suggest that the lightcurves can be explained by a 
flaring behaviour of the power-law component, while the ther- 
mal component does not vary on these short timescales, and 
represents the baseline flux seen in the soft band. 

To test our hypothesis, we analysed the lightcurve in the 
0.2-12 keV band (Fig. 3, bottom): we divided the exposure 
into a "low" interval, consisting of all the 50-s bins in which 
the count rate was < 1 ct s _1 , and a "high" interval, for rates 
> 1 ct s" 1 . We then extracted and analysed the combined EPIC 
spectra from the two sub-intervals; as before, only the 0.3-12 
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Table 5. As in Table 3, but with sub-solar abundances (Z = 
0.07Z o ) for the intrinsic absorber. 



Parameter Value in 2001 Value in 2002 Value in 2003 



model: wabsoai X tbvarabs x bknpo 





"•-0.1 


< 0.2 


< 0.2 


r, 

E bl (keV) 

r 2 

K bp (10- 3 ) 


76 +0 16 

"•'"-0.48 

55 +0 05 
3 20 +003 
3 4 +0 - 2 


o.58 + ° : « 

"•""-0.06 

2 89 +o n 

z -°'-0.07 

2 +3 3 

■^•"-0.6 


1 iQ+0.68 
55+ 05 

9 Q9+0.06 
2 -°'-0.06 

6 2 +0 - 6 

"•^-0.2 


xl 


1.78(244.0/137) 


1.09(99.2/91) 


0.97 (83.2/86) 


/0.3-12 (t) 


t 9+0.2 


9 -2+0.1 

z..o_ [ 


2 +01 

^■"-0.1 


£fl.3-12 <C) 


0.9 


0.6 


0.6 


model: wabsoai X tbvarabs x (bb + po) 


kT bb (keV) 

K bb (io- 5 ) 
r 

K po (10- 4 ) 


4 +0 - 2 

"•^-0.2 

15 +001 

"■ iJ -0.02 

9 7C+0.II 
z - ' J -0.10 

5 3 +as 


1 +0 1 
13 +001 
1 3 +0 ' 2 

-0.2 

2 6r - 06 

5 +0 5 

J - -0.2 


1 4+0.7 

12 +002 
"•^-0.02 

1 6 +L7 
2 67 +ai8 

■ t,o/ -0.17 

4 8 +L0 


xl 


1.12(153.6/137) 


0.93(84.3/91) 


0.91 (78.3/86) 


/0.3-12 <b) 


T Q+0.2 


Q Q+0.2 

-'•-'-0.1 


T Q+0.9 

J -0.4 


ifl.3-12 <C) 


1.1 


0.9 


0.9 



a in units of 10 21 cm" 2 

b in units of 10~ 12 erg cm -2 s 

c in units of 10 40 erg s _1 



keV band was used for spectral fitting 1 . This choice of thresh- 
old is purely arbitrary: we needed to have enough counts in 
both "states" for a meaningful spectral analysis. We examined 
whether the change in flux between the peaks and the troughs 
could simply be modelled with a constant factor, or whether 
there was also a significant spectral change. We find (Fig. 6; 
Tables 7, 8) that the blackbody component is consistent with 
being constant between the two states: the flux variations are 
consistent with being entirely due to a change in the power- 
law normalisation, in agreement with our initial hypothesis. In 
conclusion, we suggest that the ULX exhibits soft dips over 
timescales of ~ 10 2 s in which the power-law component is 
suppressed. 

We also fitted the two flux-dependent spectra with a bknpo 
model (Table 9). The model provides an equally good fit for 
this particular observations, as shown earlier (Tables 3, 5, 6). 
However, its physical interpretation is less straightforward. It 
requires a very high power-law index below the break (;Ti ss - 1 

1 The pn and MOS spectra were restricted to the same good-time- 
interval, for consistency with the lightcurve analysis. However, the 
exposure time is slightly longer in the MOS cameras (3.3 ks versus 
3.0 ks for the pn, Table 1) because of the higher live-time fraction 
(a 99% of the good-time-interval for the MOS CCDs, as opposed to 
a 91% for the pn CCDs). 



Table 6. Comparison between the absorbed broken power-law 
model fits ("model 1") and the absorbed blackbody plus power 
law fits ("model 2"), over the 5 observations. and are the 
null hypothesis probabilities for the two models. p2 is the prob- 
ability that model 2 is significantly better than model 1, from 
the Williams & Kloot (1953) test. (For the 2003 observation, 
neither model is significantly better than the other.) 



Fit stats 


31/07/01 


08/08/01 


24/08/01 


29/07/02 


28/01/03 






Z = 








X v , 1 


1.16 


1.45 


0.99 


1.09 


0.97 


X v ,2 


1.00 


1.09 


0.86 


0.98 


0.99 


pN 
r 1 


0.09 


1.9 x 10^ 4 


0.50 


0.26 


0.57 


pN 


0.47 


0.21 


0.88 


0.54 


0.50 


P2 


93.1% 


> 99.99% 


96.8% 


87.5% 




Z = 0.07Z o 


X v , 1 


1.15 


1.43 


0.99 


1.09 


0.97 


xl, 2 


1.00 


1.10 


0.87 


0.93 


0.91 


pN 
' 1 


0.11 


3.5 x 10~ 4 


0.51 


0.26 


0.57 


pN 
r 2 


0.47 


0.19 


0.86 


0.68 


0.71 


Pi 


93.9% 


99.9% 


97.4% 


99.4% 






0.5 1 2 5 10 

channel energy (keV) 



Fig.l. Coadded EPIC spectra of the ULX in 2001 Jul-Aug 
(blue, with its best-fit blackbody plus powerlaw model), 2002 
July (red), 2003 Jan (green). Also overplotted: the best-fit 
model for the 2001 spectrum, in photons cirT 1 s~' keV -1 (ar- 
bitrarily shifted for clarity). 

below £bi ~ 0.55 keV), thus making the soft component prac- 
tically indistinguishable from a thermal component over the 
XMM-Newton/EPIC energy band. In this scenario, the flaring 
behaviour would correspond to a flattening of the high-energy 
power-law: from T2 ~ -3.2 in the low state to T2 ~ -2.8 in the 
high state Finally, we show (Fig. 7) that the difference between 
low and high states cannot be due simply to a normalisation 
factor. The low state is indeed softer (consistent with our inter- 
pretation). 
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wabs*tbvarabs*bknpo 

abs*tbvarabs*(bb+po) 





0.5 1 2 5 10 

channel energy (keV) 

Fig. 2. Coadded EPIC spectra of the ULX on 2001 Aug 08, 
with the best-fit bb+po and bknpo models (green and red 
curves, respectively). A metal abundance Z = 0.07 Z Q was as- 
sumed for the intrinsic absorber. See Table 4 for the fit param- 
eters. 




2001 Aug 25 




2002 Jul 29 -E 



Time (s) 

Fig. 3. Coadded EPIC lightcurves of the ULX over the four 
observations of 2001-2002. 





Fig. 4. Coadded EPIC lightcurves of the ULX in 2003 Jan 
(Rev. 574), rebinned to 50 s. Top panel: 1.5-12 keV band; 
medium panel: 0.2-1.5 keV band; bottom panel: full 0.2-12 
keV band. 




10 J 0.01 0.1 

Frequency (Hz) 

Fig. 5. Power density spectrum of the source in 2003 Jan, nor- 
malised to give the squared rms fractional variability per unit 
frequency; the expected white noise level has been subtracted. 
Only the first, uninterrupted 3 ks of the exposure have been 
used for this spectrum; however, very similar results are ob- 
tained when we used the whole exposure. 



To investigate further the characteristic variability 
timescales, we examined the power spectral density for the 
2003 lightcurve; we normalised the Fourier transform to 
be the squared rms fractional variability per unit frequency. 
Firstly, we considered in our analysis only the first unin- 
terrupted 3 ks of the exposure (Fig. 4); then we considered 
the power spectrum for the whole exposure, filling the two 
data gaps with a running mean. We obtained similar results. 
The power spectrum for the full-band lightcurve (Fig. 5) 
has a slope a = -1.3 ± 0.2 over two decades in frequency 
(3 x 10~ 3 < v $ 0.3 Hz); it peaks at v b * 2.5 x 10~ 3 Hz, and 
is flat (consistent with a = 0) at lower frequencies. The power 
spectrum for the hard-band lightcurve has similar features, 
with a break at Vb ~ 3 x 10~ 3 Hz and a slope a = -0.9 + 0.4 at 
higher frequencies. 



4. Conclusions 

Our spectral analysis, based on five XMM-NewtonfEPIC obser- 
vations between 2001 July and 2003 January, shows a thermal 
component with kT\,b x 0.12-0.14 keV, significantly detected 
at all epochs, in addition to a soft power-law (F as 2.6-2.9). A 
broken power-law model, as suggested in Kaaret et al. (2003), 
is acceptable at some epochs, but is not consistent with the full 
series of observations, and is generally worse than a blackbody 
plus power-law model or a Comptonised blackbody model. 

Ruling out a simple broken power-law model does not 
rule out the possibility that the X-ray emission is due to 
inverse-Compton emission from a jet, enhanced by relativis- 
tic Doppler boosting (Kaaret et al. 2003). It was suggested 
in Georganopoulos (2002), and Georganopoulos & Kazanas 
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Table 7. Best-fit parameters for a bb+po simultaneous fit to the 
combined EPIC spectra in the "low" and "high" intervals of the 
2003 Jan observation. As before, the quoted errors are the 90% 
confidence limit and «H,Gai = 5.7 x 10 2() cirT 2 . We assumed a 
solar abundance for the intrinsic absorber (Z = Z ). 



Parameter "Low" value "High" value 



model: wabsoai X wabs x (bb + po) 



n„ (10 21 cm- 2 ) 




(1.4!»|) 


(1.4t°J) (a) 


kT bb (keV) 






12 +001 


Ab b (10- 5 ) 




3.7^ 


~> Q+4.8 


r 




2.86!°; 42 


2.68!°;!? 


K po (ltr 4 ) 






6.4!|; 2 


xl 




0.76(137.8/182) 


/o.3-i2,obs (10~ 12 erg cm -2 




i r - 2 

1 ' 1 -0.4 


2 +0 - 2 

^'"-0.4 


/o.3-i2(10- 12 ergcm- 2 s" 


) 


3 6 +(u 

—1.9 


5 5 +03 

-2.3 


/o b 3-i2(10~ 12 ergcm- 2 s" 


) 


2 r°- 7 


2 3 +L1 


/ p °_ 12 (10- 12 ergcm- 2 s" 


) 


1 5 +L2 


3 1 +L3 


Lo.3-12 (10 40 erg s" 1 ) 




1.0 


1.8 



a assumed to be equal for the two states 



Table 8. As in Table 7, but with Z = 0.07 Z c 



Parameter "Low" value "High" value 



model: wabsoai X tbvarabs x (bb + po) 



n H (10 21 citT 2 ) 




(1-611) 


(1.6tl) (a) 


kT bb (keV) 




o-iiS 




Abb (10- 5 ) 




1 9 +1 8 


1 7 +L9 

'-■ '-0.8 


r 






2 66 +018 

^• ou -0 .19 


A- po (10- 4 ) 




3 0+°- 9 

j.u _ 07 




xl 




0.72(130.9/182) 


/o.3-i2,ob s (10 -12 erg cm -2 




1 l +0 - 2 


2 +0 - 2 


/o.3-i2(10- 12 ergcm- 2 s" 


) 


2 6 +0 - 8 

z.u 03 


4 +0 - 9 


/o b 3-i2(10"' 2 ergcm- 2 s" 


) 


1 l +0 - 5 


1 o +0 - 7 

i.o_ j 


/ p °_ 12 (10- 12 ergcm- 2 s" 


) 


1 5 +03 


3 +03 

J - -0.3 


Lo.3-12 (10 40 erg s" 1 ) 




0.7 


1.1 



a assumed to be equal for the two states 



(2003), that a more realistic jet spectrum would have a cur- 
vature around the break energy, thus mimicking a soft thermal 
component. Spectral analysis over a larger energy range is re- 
quired to test this possibility. 

In any case, the X-ray spectrum is consistent with a more 
massive accreting BH in a high/soft or very high state. Given 
its inferred isotropic luminosity of a 10 40 erg s -1 in the 0.3- 
12 keV band, a mass 100M o is required to satisfy the 
Eddington limit. Thus, the ULX in NGC 5408 appears very 
similar to a group of "canonical" ULXs in nearby dwarf and 




0.5 1 2 

channel energy (keV) 



Fig. 6. Coadded EPIC spectra of the ULX on 2001 Jan 28, sep- 
arated into the "low-flux" and "high-flux" states (blue and red 
datapoints, respectively). See text for details. The spectra are 
well fitted by a bb+po model: we argue that the blackbody com- 
ponent is consistent with being unchanged in the two states, 
while the variability is due to the power-law component. A 
metal abundance Z = 0.07 Z Q was assumed for the intrinsic 
absorber. See Table 8 for the fit parameters. 




0.5 1 2 5 10 

chonnel energy (keV) 



Fig. 7. Same datapoints as in Fig. 6. The fitted bb+po model 
shows that the change in flux between low and high states can- 
not be due simply to a normalization factor: the spectrum in the 
low state is softer. 



spiral galaxies, with a low-temperature thermal component at 
kT bb ~ 0.1 keV and emitted luminosity ~ (1-2) x 10 40 erg 
s~'. In all these systems, if the soft thermal component is in- 
terpreted as the emission from the inner part of a standard ac- 
cretion disk (Shakura & Sunyaev 1973), its temperature would 
imply masses > 1O 5 M , unphysically high for a non-nuclear 
BH. However, it was also suggested (King 2003) that these 
ULXs may be accreting BHs during a super- Eddington ac- 
cretion phase: they could have a luminosity close or even a 
factor of a few above the Eddington limit, accompanied by 
a radiatively-driven, Compton-thick outflow from the accre- 
tion disk. In this scenario, the soft thermal component would 
come from the photosphere of the outflow. Hence, the BH mass 
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Table 9. As in Table 6 and 7, for a bknpo model. We assumed 
Z = O.O7Z . 



Parameter "Low" value "High" value 



model: wabsoai X tbvarabs x bknpo 



n H (10 21 cm- 2 ) 


(< 0.6) 


(< 0.6) (a) 


r, 


-1 1 +24 

1 - 1 -1.8 


-1 3 +2 -° 
±.j_ 02 


£ br (keV) 


55+ 009 

-0.05 


56+ 006 


r 2 


o nc+0.21 
- , ' iJ -0.16 


9 7Q+0.12 
' 7 -0.08 


K bp (10- 3 ) 


A Q+10.1 


7 1 +13 - 5 
'•-5.0 



Xl 0.75(137.1/182) 

Zo.3-i2.obs (10- 12 erg cm- 2 s" 1 ) 1.0*°;] 2.0+°;] 

/0.3-12 (10- 12 erg cm- 2 s" 1 ) 1.4+"; 2 2.4+°;] 

io.3-12 (lO^ergs- 1 ) 0.4 0.7 



a assumed to be equal for the two states 



would not need to be higher than ~ 2O-5OM , still consistent 
with a stellar origin. 

We note that this ULX is located in a very metal-poor 
dwarf galaxy (Z « O.O7Z ). An association between ULXs and 
metal-poor environments was suggested by Pakull & Mirioni 
(2002). Low metal abundance implies a reduced mass-loss 
rates in the radiatively-driven wind from the O-star progeni- 
tor (M w ~ Z 85 : Vink et al. 2001; see also Bouret et al. 2003); 
this leads to a more massive stellar core, which may then col- 
lapse into a more massive BH, via normal stellar evolution. 
Hence, low abundances may explain the formation of isolated 
BHs with masses up to ss 5OM (ie, up to « half of the mass of 
the progenitor star). 

We found a moderate long-term spectral variability be- 
tween the various epochs of the XMM-Newton observations, 
with a ratio between the maximum and minimum fluxes of 
« 1 .4. This is consistent with the long-term behaviour of the 
source from previous Einstein, ROSAT, ASCA and Chandra ob- 
servations (Kaaret et al. 2003). Our preliminary timing analysis 
shows rapid variability and flaring-like behaviour, particularly 
in the 2003 Jan observation. The variability is larger in the hard 
band (E > 1.5 ke V), where the flux is consistent with zero dur- 
ing the dips. 

The bb+po spectral model offers a clean, simple interpre- 
tation of the flaring behaviour. In this scenario, we have inter- 
preted the lightcurves as the combination of a rapidly variable 
power-law component, plus an additional, non-variable ther- 
mal component below 1.5 keV, consistent with our spectral 
analysis. The power-law flux appears to be suppressed during 
the dips, when the X-ray spectrum becomes softer; this is sup- 
ported by flux-dependent spectral analysis. The power-law X- 
ray emission in accreting BH binaries is generally explained as 
inverse-Compton scattering of soft disk photons in a hot corona 
or Compton cloud, located either above the inner part of the ac- 
cretion disk, or as a quasi-spherical region inside the truncation 
radius of the disk. 



A flaring behaviour with soft dips or transitions has 
also been observed in Galactic microquasars such as GRS 
1915+105 (Naik et al. 2001) and, in one case, XTE J 1550-564 
(Rodriguez et al. 2003). It was interpreted (Vadawale et 
al. 2003) as evidence that the matter responsible for the 
Comptonised component is recurrently ejected from the in- 
ner regions; the ejected matter may be responsible for the 
optically-thin synchrotron radio emission in those systems. A 
more detailed discussion of this hypothesis, and of whether it 
may explain the steep-spectrum radio emission detected by the 
Australia Telescope Compact Array at the position of the ULX 
(Kaaret et al. 2003; Stevens et al. 2002), will be addressed in 
further work. 

The power density spectrum for this ULX is flat below 
« 2.5 mHz, and has a slope of ^ -1 at higher frequencies. 
This is similar to what is found in "canonical" Galactic BH 
candidates and AGN: it is believed that the break frequency 
is inversely proportional to the mass of the accreting BH (eg, 
Belloni & Hasinger 1990; Nowaket al. 1999; Uttley et al. 2002; 
Markowitz et al. 2003; Cropper et al. 2004). For example, 
the break in the power density spectrum of Cyg X-l occurs 
at v ~ 0.4-0.04 Hz in different accretion states (Nowak et 
al. 1999). This suggests that the mass of the BH powering the 
ULX in NGC 5408 is ~ 1O 2 M . This is also in good agree- 
ment with the inferred luminosity, if we require that it does not 
exceed the Eddington limit. Such a high mass would rule out 
(cf. Fig. 3 of Kaaret et al. 2004) an association of this ULX with 
the young star clusters in the starburst region of NGC 5408, lo- 
cated « 300 pc to the north-west. However, we are aware that 
the relation between break frequency and BH mass is very un- 
certain, and the power-density spectra are often more compli- 
cated than a simple broken power law (see, eg, the case of GRS 
1915+105: Morgan etal. 1997). 

In conclusion, we argue that both the spectral and timing 
results from our XMM-Newton study are consistent with the 
possibility that this ULX is at the upper end of the mass range 
for BHs of stellar origin, a few times more massive than "stan- 
dard" BH candidates such as Cyg X-l. In that case, its luminos- 
ity would be very close or, more likely, slightly higher than its 
Eddington limit, and a strong outflow would be expected from 
its accretion disk. Doppler boosting in a relativistic jet is not re- 
quired to explain the observed properties, although we cannot 
rule it out. 
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